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1. Introduction 


The  gravity  potential  of  the  earth  (W)  can  be  decomposed  into  the 
gravitational  potential  (V)  and  the  rotational  potential  ($): 

W = V + $ 


As  V is  harmonic  outside  the  surface  of  the  earth,  it  can  be  expanded  into  a series 
of  spherical  harmonics  in  this  region  (we  assume  that  the  series  is  convergent): 

00  n 

(1*1)  V=  {1  + V Y <c  n»  cosmA+ ?nB  sinmA)  PnB(sin  <p)  j. 

n=3  n=0 

where 

G Newton's  constant  of  gravitation 


M mass  of  the  earth 

( r, <p , A ) geocentric  spherical  coordinates  of  the  computation  point 
Pn»()  normalized  associated  Legendre  polynomial 


C n»  » S »■  fully  normalized  spherical  harmonic  potential  coefficients 
a equatorial  radius  of  the  defined  earth  ellipsoid 

In  the  study  of  the  geopotential  field  it  is  most  convenient  to  subtract  a selected 
reference  field: 

2 

u = ~~  i1  + ca°  (~)  Pao  (sm<P)  + c40  (-|-)  P40(sin<p)j-  + § 


A A 

where  Cao  and  C4o  are  selected  coefficients,  from  the  potential  W.  The  residual 
potential  (T)  is  the  so  called  disturbing  potential; 

00  a 

(1*2)  T = W - U = ^ (rj  / (Cn*  cosm^  + sn«  sinmA)PBB  (sin <p) 

n=  a B = o 


"i  9 01  12  004 


where 


Cao-  Cac-Cao 


C i o C /|  o C A ( 


c„.  = CnB,SnB=SnB  for  (n,m)  ^ (2,0)  or  (4,0) 


Now,  by  inserting  (1.2)  into  the  spherical  approximation  of  "the  fundamental 
boundary  condition"  (He iskanen  and  Moritz,  1967,  p.  88): 


Ag  = . |T  . IT 
or  r 


we  obtain  the  gravity  anomaly  (Ag)  in  terms  of  the  potential  coefficients: 

oo  n 8 

(1.3)  Ag  = y 5>1)£  (Cnl  cosmA  + S„,  sinmA)^)  P„,  (sin<0) 


n=2  o— 0 


where 


y = GM/az 


It  is  obvious  from  formula  (1.1)  that  the_coefficients  C nB  and  S’,,,  are  dependent 
on  the  choice  of  a,  but  that  C na  a"  and  S nB  an  are  invariant  quantities.  On  the 
other  hand,  for  a fixed  y the  gravity  anomaly  in  (1.  3)  is  independent  of  the  choice 
of  a,  whenever: 

CnBan+a  and  SnB  an+a 


are  invariant.  This  means  that  for  two  different  a-values  (ax,  aa)  and  y = constant 
the  gravity  anomalies  in  (1.  3)  are  the  same,  if  (with  obvious  notations): 

r C / O -\B  + a f C nil 


= (M  {.} 

lsB.Jn  W ^SnlJ.a 


This  relation  will  be  useful  in  the  following  application  of  estimating  the  coefficients 
C n i and  §"nB . 
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2.  Covariance  Functions 


Suppose  that  the  gravity  field  is  harmonic,  homogeneous  and  isotropic. 
Then  the  spatial  covariance  function  of  the  point  free-air  gravity  anomalies  Agi 
and  Agj  is  given  by  (Moritz,  1972,  p.  89): 


(2- !)  c u = cov  (Ag  j , Ag  j)  = Y c „ sn+a  pn  (cos  0 u ) 

n=a 


where 


s = re2/r  j rj 

re  = radius  of  the  Bjerhammar  sphere 

rk  = r,  + hk,  k = i,  j 

r„  = radius  of  mean  sea  level 

h j , hj  = elevations  of  points  P,  and  Pj 

cn  = degree  variances  of  Ag 


The  corresponding  covariance  function  of  the  mean  gravity  anomalies  Ag  j and 
Ag  j , where: 


**'m  sk  If  ABd° ! k“'i 

act„ 


is  given  by 


(2.2)  Jll-COT(Ag„AgJ).siaj  rf  JJ  C.dOdO 

Aa(  AOj 

In  the  same  way  we  obtain  the  following  cross-covariance  function  between  Ag( 
and  Agj: 


(2-3)  ClJ.  _i.  JJ  ClJ 

40, 

In  these  formulae  A a is  a part  of  the  unit  sphere  a.  For  small  regions  Adj  and 
Affj  we  may,  without  loss  of  significance,  assume  that  r,  and  r3  are  constants. 
Then  (2.2)  and  (2.3)  become; 

c- s'*“  35^7  II JJ  p.«*»*„)dodo 

»=»  A a,  Affj 


(2.2  ) 


The  common  way  to  determine  the  spherical  harmonic  coefficients  of  the 
earth's  gravity  field  from  terrestrial  data  is  by  means  of  integration  of  mean 
gravity  anomalies  over  a mean  earth  sphere.  See  for  example  Rapp  (1977a).  Due 
to  the  orthogonal  functions  the  coefficients  are  obtained  directly  in  the  integrations. 
Rigorously,  it  is  required  that  the  gravity  anomalies  are  located  on  the  sphere  of 
integration,  but  in  reality,  due  to  the  variation  of  the  elevation  of  the  terrain,  this 
is  not  the  case.  This  terrain  deviation  can  be  corrected  for  by  adding  the  Molo- 
denshy  Gi  term  to  the  mean  anomalies.  Pellinen  (1962)  has  indicated  that  the 
neglect  of  this  term  can  cause  errors  in  the  low  degree  coefficients  of  10  to  20 
percent.  Numerical  results  of  Rapp  (1977a)  agree  with  this  error  estimate,  but 
the  results  were  based  on  a number  of  assumptions  relating  the  Gi  term  to  ihe 
terrain  correction. 


In  practice,  the  computation  of  the  Molodenshy  correct.,  term  for  the  terrain 
may  be  very  difficult  and  laborious,  especially  in  areas  with  rapidly  varying  topog- 
raphy. It  is  therefore  of  interest  to  find  a technique  for  the  determination  of  the 
potential  coefficients  that  does  not  include  the  computation  of  the  Gi  terms,  yet 
retains  the  rigor  of  that  procedure.  One  such  method  was  given  by  Rapp  (1977b), 
where  least-squares  collocation  was  used  for  an  upward  continuation  of  the  terres- 
trial mean  anomalies  to  a bounding  sphere.  Once  the  anomalies  are  given  at  the 
sphere  the  integration  can  be  applied  strictly  for  the  determination  of  the  potential 
coefficients.  Rapp  (1977b)  found  that  the  neglect  of  Gx  caused  errors  less  than 
7.  5%  for  harmonics  up  to  degree  40. 


In  the  present  study  the  idea  is  to  estimate  the  potential  coefficients  directly  by 
applying  least-squares  collocation  to  the  terrestrial  mean  gravity  anomalies.  The 
integration  is  then  taken  care  of  in  the  cross -covariance  matrix.  The  advantage  of 
using  such  a method  is  that  the  terrain  correction  is  easily  included  and  that  the 
various  accuracies  of  the  mean  anomalies  can  be  taken  into  account,  which  is  not 
obvious  in  the  integration  approach. 


In  collocation  a physical  quantity  V may  be  predicted  from  a vector  of  (mean) 
gravity  anomalies  _Ag  by  the  relation: 


(1.5) 


V = cyT  (G  + D)_i  Ag 


where 


The  prediction  errors 


CyT  = cross -covariance  matrix  (V,  Ag) 
C = auto-covariance  matrix  ( Ag,  Ag) 

D = error  covariance  matrix 
are  estimated  by 


(1.6) 


m»  — Co  - CyT  (C  + D)  1 Cy 


Where  Co  is  the  variance  of  V prior  to  prediction.  For  further  details  on  these 
basic  formulae  see  Moritz  (1972). 


The  collocation  technique  requires  that  the  relevant  covariance  functions  are 
known.  In  the  next  section  we  are  going  to  study  these  functions  for  the  present 
application. 


-4- 


and 


1 


(2.3  ) 


CO 


[ Pn  (COS  (|>,j)  do 
AOj 


Formulae  (2.2  ) and  (2.  3 ) are  very  laborious  to  compute  in  practice  due  to  the  many 
numerical  integrations.  Approximate  mean  covariance  functions  may  be  obtained  by 
using  the  /?„  function  of  Meissl  (1971,  p.  23): 


1 

1-cos  i/)0 


1 

2n+l 


[ P n— 1 (COS  ii0)  " ^n  + 


(cos  ib0)] 


where  ib0  is  the  radius  of  a circular  cap  of  area  equal  to  the  relevant  blo-k  size  of 
the  mean  anomalies.  Then,  the  above  covariance  functions  become  (approximately) 
(cf.  Figure  2 a-b): 


(2*2  ) C jj  = 7 C „ s"+S  Pn  (COS  0 t j ) 

n = 2 


and 

00 

(2. 3 ) c,j=V  cnsn+a£n  P„(cos0u) 

n=  a 

As 


£n  - 0 as  n-00 


it  is  usually  sufficient  to  truncate  the  series  (2.2'")  and  (2.3")  at  a few  hundred  degrees 
(dependent  upon  the  block  size)  without  loss  of  significance. 

We  also  give  the  autocovariance  relations  between  the  potential  coefficients. 

In  Moritz  (1970)  the  relations  are  given  for  the  anomaly  coefficients  (ani,  bnB).  As 
anl  and  bnB  are  related  to  the  potential  coefficients  C n,  and  Sn.  according  to: 


} 


where  y is  the  mean  gravity  at  sea-level,  we  obtain  from  Moritz  (ibid.): 
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cov  (C  na , Sp  q)  = 0 in  any  case. 


The  above  covariance  relations  will  be  useful  in  our  application  . collocation. 

Finally,  we  like  to  mention  that  the  mean  covariance  functions  (2.2  ) and 
(2.  3 ) can  in  addition,  be  approximated  by  the  corresponding  point  covariance  functions 
simply  by  increasing  the  radii  r4  and  rj  by  a feasible  constant.  This  type  of 
smoothed  covariance  functions  was  discussed  in  Tscherning  and  Rapp  (1974,  Section  10) 
and  Schwarz  (1976,  Section  7). 


3.  Application  of  Collocation 

We  assume  that  the  external  gravity  field  of  the  earth  may  be  expanded  into 
a series  of  spherical  harmonics  at  a sphere  of  radius  R.  Then  we  have  [cf.  (1.3)]: 


00  n+  a n 

la)  Agn  = y V (n-  1)  ^ V (C  n,  cosmA  + SnB  sinm  A)  Pn„  (sin  <p) 


where 

y = GM/a2 

re  = radius  of  the  Bjerhammar  sphere 


(3.  lb) 


{.}  - 

LS..  J r n 


» JrB  4ny(n-l)(rB/R)  o 


n cosmA. 

Ago  Pr..  (sin  <p)-|  da 

w s in  m A J 


Formula  (3.  lb)  is  the  basic  equation  we  are  going  to  use  for  estimating  the  potential 
coefficients.  The  coefficients  determined  by  (3.  lb)  are  independent  of  the  choice  of 
R.  Thus  by  choosing  R as  the  radius  of  the  Brioullin  sphere  (bounding  all  the  mass 
of  the  earth)  we  have  a theoretically  most  attractive  situation,  because  the  series 
expansion  (3.1a)  is  rigorously  convergent  at  this  sphere  (cf.  Sjoberg,  1977).  The 
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standard  representation  of  the  coefficients  C „ and  S ni  in  the  literature  is  with 
reference  to  the  equatorial  radius  of  the  earth  (a)  (and  not  to  the  Bjerhammar 
sphere  as  in  (3.  lb)).  The  conversion  from  (3.  lb)  to  this  representation  was 
given  in  (1. 4): 


n + 3 c 1 


(3. 2) 


{,!.-(?)  V) 

am  d ^ n» 


rB 


Theoretically,  the  point  gravity  anomaly  AgR  in  (3.  la)  can  be  estimated 
from  a vector  of  mean  gravity  anomalies  by  means  of  formula  (1.3): 


(3.3) 

where 


Ags  = crT  (C  + Df1  Ag 


cr  = cross-covariance  matrix  (Ag , Ag) 


The  element  (csT),j  is  given  by  (2.3')  with  r ,=  R and  C tJ  is  given  by  (2.2').  The 
error  covariance  matrix  D can  be  estimated  bv  the  diagonal  matrix  formed  by  the 
a priori  estimated  mean  anomaly  variances. 

By  inserting  (3.  3)  into  (3.  lb)  we  arrive  at: 


r c nm  ■> 

(Cc) 

i ^ r 

1 T f 

^ nB 

'‘Cs  J 

(3.4) 

where  the  elements  of  cc  and  cs  become; 

(3.5) 


r( c c ) U •> 

1 

» , _ cosmA. 

HCslijJ 

4TTy(n-l)(re/R),,+a  1 

JJ  Cl'p"<sln,0){Si»mx} 

rr 

da 


We  are  going  to  simplify  (3.5)  by  taking  into  account  the  orthogonality  of  the  spherical 
harmonics  (Heiskanen  and  Moritz,  1967,  p.  29).  Using  the  notations; 


f Rn»  (V.M  _ r cos  m A.-V 

{_  f = Pn.  (Sin0){  } 

Ls„.  (<p,A)j  oi~  ™ ' 


j-  cos  m A> 
•sin  m A' 


we  obtain  from  (2.3  ); 
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OD 

(3.6)  4^  If  ?,J  Sn.«0,X)dCT  = ~ Jj  Y CnSB+S  J j Pn  (COS  0 „)  Rot  (0 , X)  do  do 

J AOj  »=a 


sn 

2n+l 


“ ^ JJ  d0 


A completely  analogous  result  is  obtained  for  SnB.  Inserting  the  result  of  (0.6)  into 
(3.5)  we  finally  obtain: 


r(cc)u,  c„  l pf-  C””  ‘"'"1  j 

(c s)ijl  ~ y(2n+l)(n-l)  (rj)  Act3  Jj  1 aa  ’in0)  (sin  mX  J 


cos  mA, 


where 


. cos  m X 


j'J  da 


^ N 

= — f Pn»  (sin<P)  cos 

sin0N-sin0s  *) 


r 1 if 

I (sinr 
<p  d 0 x / 

I (cos  1 


1 if  cos  m X with  m = 0 
(sinmAt-sinm  X*)/m  if  cos  mX 
with  m ^ 0 

(cos  m A.w  - cos  m Xt)/m  if  sin  m X 
with  ra^O 
0 if  sinmX  with  m = 0 


0s , 0N  , Xe,  Xw  = geocentric  latitudes  and  longitudes  of  the  corners  of 
the  block  A 0 s 


If  c jj  is  approximated  by  (2.  3~)  the  formula  analogous  to  (3.  7)  becomes: 

f(c‘>u,  cn0n  /re\n+s  _ rcosml'i 

(3‘?  ] MCs)ij}  y(2n+l)(n-l)  (rt)  <P)  Vsin  mX  ^ 


where 


0 = (0N  + 0s)/2 

X = (Xe  + X*)/2 


\«r< 


4.  Error  Analysis 

Suppose  that  the  gravity  anomaly  Agi  used  in  formula  (3.1b)  has  the  error 
( i . The  error  propagation  to  C M is  accordingly: 


dt  n«-  IT  f|  Rn,  ((0|X)  d a 

4ny(n-l)(re/R)n+s 


and  the  mean  square  error  of  C na  becomes: 


me*.  = M [dC„,}  = 


^y(n-l)(rB/R)n*a}  M 1 ^ Rn»(0i  .A,)  Rn.  (Vj.Aj)}  } 


where 


Ou  = M {e,Cj}  = ^ JJ  e,  Cj  d a 

l/i  i j = const. 


As  Ag  of  formula  (3.  lb)  is  estimated  by  means  of  collocation,  au  is  the  prediction 
covariance  of  that  method.  Moritz  (1972),  formula  (3-39)  with  A = 0 gives 


<*!j  = C ,j  - c,T  (C  + D)-i  cj 


where  the  covariances  of  the  right  member  are  those  defined  in  (3.3).  Inserting  (4.2) 
into  (4.1)  and  carrying  out  the  integrations,  we  finally  obtain: 


(4.3a) 


— 2 C 

mcB.  = 11 r— - - ccT  (C+  D)~l  cc 

(2n+l)(n-l)ay 


where  the  elements  of  cc  are  given  by  (3.  7).  In  the  same  way  we  obtain  the  mean 
square  error  for  S„,  : 
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(4.  3b) 


ms„. 


c 


* 


n 

(2n+l)(n-l)a  ya 


csT  (C  + D)'1  cs 


We  notice  that  (4.3a-b)  give  exactly  the  error  estimates  we  would  expect  in  collocation 
with  the  prediction  formula  (3.4).  Cf.  formulae  (1.5)  and  (1.6). 

5.  Computations 

A global  coverage  of  416  10°  equal  area  free  air  gravity  anomalies  were 
available  as  input  data.  (These  anomalies  had  been  derived  using  the  data  of  Raop 
(1977a)  and  the  methods  described  by  Hajela  (1975).)  In  all  commit  I covariance 
functions  the  degree  variances  implied  by  the  subroutine  COVA  t Tseheming  and 
Rapp  (1974)  were  used  with  ca  = 7.5  mgal2.  The  cross-co\  riance  functions  were 
computed  by  numerical  integration  according  to  formula  (3.7)  [derived  from  (2.3)j. 
However,  for  the  auto-covariance  function  we  felt  that  it  was  unreasonable  to  use  the 
corresponding,  very  laborious  formula  ( 2.2 ').  Instead,  we  tried  two  different 
approximate  formulae.  In  order  to  save  computer  time  the  auto-covariance  function 
for  the  mean  anomalies  for  zero-elevations  (h,=  hj  = 0)  was  stored  in  a table  and 
the  current  values  were  interpolated  among  the  tabulated  values.  When  elevation 
information  was  included  in  the  process  the  table  was  used  only  to  determine  the 
auto-covariances  between  ocean-block  mean  anomalies. 


The  following  reference  data  were  anticipated: 


a = 6378140  m 

Cao=  -484. 198  x 10*6 

C 4 o = 0.790333  x 10“6 

f = 1/298.247 


In  a first  test  of  the  prediction  formulae  (3.4)  and  (4.3  a-b),  the  auto- 
covariance function  was  approximated  by  the  corresponding  point  covariance  function 
at  the  best  fitting  elevation  (h  , = hj  = 142.  3 km,  cf.  the  end  of  section  2).  Figure  1 
indicates  a fair  agreement  between  this  covariance  function  and  the  one  implied  by  a 
numerical  integration  of  the  point  covariance  function  over  10°x  10°  blocks  (around 
the  equator).  However,  the  prediction  result  was  poor  and  especially  the  error 
estimates  were  useless  (because  of  negative  variances !).  It  was  apparent  from  the 
test  that  choice  of  mean  covariance  function  was  very  critical,  especially  for  the 
outcome  of  the  prediction  errors. 

Second,  the  auto-covariance  function  was  computed  according  to  formula  (2.3  ). 
It  was  found  that  the  series  could  be  accepted  when  truncated  at  degree  200  (except  for 
= 0,  where  n = 200  is  sufficient).  The  very  good  agreement  between  the  series 
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(2. 2") and  the  numerically  integrated  10°x  10°  covariance  function  [formula  (2.2')] 
is  illustrated  in  Figure  2a-b.  In  Figure  2a  we  can  distinguish  between  the  two 
functions,  in  Figure  2b  they  coincide. 

The  spherical  harmonic  coefficients  were  determined  to  degree  20,  both 
with  mean  elevations  set  to  0 (Table  A.  1)  and  with  the  inclusion  of  approximate  mean 
elevations  (Table  A.  2).  (The  approximate  elevations  were  computed  from  the  5° 
anomaly  data  reported  in  Rapp  (1977a).)  We  also  computed  the  coefficients  b> 
collocation  for  zero  elevation  and  no  error  covariance  matrix  included  (Table  A.  3). 

The  tabulated  coefficients  refer  to  the  sphere  bounding  the  earth  ellipsoid  with 
a = 6378140  meters  (cf.  formula  (3.2)).  For  comparison  we  determined  also  the 
coefficients  from  the  10°  mean  anomalies  by  the  integration  method  described  in 
Rapp  (ibid.)  [see  also  formulae  (6.5a-b)].  See  Table  A.  4. 

In  Table  1 we  give  the  differences  between  differeni  sets  of  potential 
coefficients.  The  relative  mean  differences  give  an  over  all  view  of  the  agreement 
between  various  sets  of  coefficients.  It  is  obvious  that  the  coefficients  implied  by 
the  5°  anomalies  (#251,  Rapp,  ibid.)  agree  slightly  better  with  GEM  9 (Lerch  et  al. , 
1977)  than  those  obtained  from  10°  anomalies.  We  also  notice  that  the  differences 
between  the  GEM  9 solutions  and  the  solutions  by  collocation  are  less  than  the 
differences  GEM  9 - Integrated  10°  anomaly  coefficients. 

In  Table  2 the  potential  coefficient  solutions  from  10°  mean  anomalies  are 
compared  relative  to  Kaula's  rule.  For  low  degree  coefficients  the  differences 
between  collocation  (1  or  2)  and  integration  is  within  6%  (except  for  n = 3).  For 
higher  degrees  the  discrepancies  increase.  In  general,  the  inclusion  of  elevation 
information  in  collocation  (Coll  1)  seems  to  change  the  estimates  a few  percent. 
However,  note  the  large  discrepancies  (10%)  for  n = 3.  When  the  noise  covariance 
matrix  is  excluded  in  collocation  (Coll  3)  the  solution  is  clearly  impared  when 
compared  with  GEM  9.  One  could  expect  that  the  collocation  3 solution  should  agree 
better  with  the  purely  integrated  coefficients  than  the  other  collocation  solutions.  Wo 
have  no  explanation  of  why  this  is  not  the  case  in  the  computations. 

In  Table  3 the  RMS  accuracy  estimates  by  degree  of  the  potential  coefficients 
are  compared.  We  notice  that  the  accuracy  estimates  by  collocation  attenuate  much 
slower  than  those  implied  by  the  integration  method.  Comparing  the  accuracy  estimates 
of  pure  integration  of  5°  and  10°  anomalies,  we  notice  that  the  10°  estimates  are 
smaller.  As  this  result  cannot  agree  with  reality,  it  is  likely  to  believe  that  the  model 
used  for  computing  the  error  estimates  is  too  rough.  The  most  obvious  reason  is 
probably  that  the  covariances  between  the  mean  anomaly  errors  are  disregarded. 
Although  this  is  the  case  also  in  the  collocation  error  estimates,  the  negative  effect 
is  less  pronounced  for  these. 

Finally,  we  compare  the  computation  time  for  the  determination  of  the  potential 
coefficients  from  10°  mean  anomalies  (Table  4).  The  integration  method  is  at  least 
30  times  faster  than  collocation. 
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Figure  2a -b.  Mean  auto-covariance  functions  for  Ag  from  subroutine 
COVA  with  ca  = 7.5  mgal2.  In  Figure  2b  the  small  scale  does  not  allow 
for  separating  the  two  covariance  functions.  In  the  larger  scale  of 
Figure  2a  the  differences  can  be  distinguished. 
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JJ 


Coll  1 = Collocation  with  Elevationa,  Coll  2 = Collocation  without  Flevatlone,  Coll  2 - Collocation  without  f leva  i Inna  »nH  without  Nolae 

Covariance  Matrix 


Table  2 


Ratios  between  RMS  Potential  Coefficient  Differences  by  Degree  Implied  by 
10°  Anomalies  and  (Cn„,Sni)  = 10-5/ns. 

Units:  Percent 


Degree 

Int- 
Coll  1 

Int- 

Coll  2 

Int- 
Coll  3 

Coll  1- 
Coll  2 

Degree 

Int- 
Coll  1 

Int- 
Coll  2 

Int- 

Coll  3 

Coll  1- 
Coll  2 

2 

3 

4 

6 

2 

12 

5 

4 

7 

2 

3 

10 

10 

9 

10 

13 

7 

7 

5 

2 

4 

5 

6 

12 

4 

14 

7 

6 

7 

2 

5 

6 

6 

5 

3 

15 

7 

7 

6 

2 

6 

4 

5 

9 

3 

16 

9 

9 

6 

2 

7 

4 

3 

4 

3 

17 

14 

14 

11 

1 

8 

6 

6 

8 

2 

18 

17 

18 

14 

1 

9 

5 

4 

5 

2 

19 

21 

22 

19 

1 

10 

4 

4 

8 

2 

20 

28 

28 

24 

1 

11 

6 

6 

6 

2 

Coll  1 = Collocation  with  Elevations 
Coll  2 = Collocation  without  Elevations 

Coll  3 = Collocation  without  Elevations  or  Noise  Covariance  Matrix 


Table  3 


RMS  Accuracy  Estimates  by  Degree 
for  Pote  *!al  Coefficients 


Degree 

#251 

Int  10° 

Coll  1 

Coll  2 

2 

1911 

1722 

2054 

1779 

3 

956 

861 

1674 

2247 

4 

637 

574 

992 

1207 

5 

478 

430 

691 

807 

6 

382 

344 

468 

438 

7 

318 

286 

443 

486 

8 

273 

245 

370 

378 

9 

239 

214 

341 

360 

10 

212 

190 

301 

300 

11 

191 

171 

285 

289 

12 

173 

155 

272 

277 

13 

159 

142 

266 

266 

14 

147 

130 

250 

248 

15 

136 

121 

243 

243 

16 

127 

113 

237 

237 

17 

119 

105 

231 

231 

18 

112 

100 

224 

224 

19 

106 

94 

217 

217 

20 

100 

89 

209 

209 

Coll  1 = Collocation  with  Elevations 
Coll  2 = Collocation  without  Elevations 
All  values  multiplied  by  lO10 
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Table  4 


Computation  Times  for  Various  Methods 


6.  An  Extended  View 

An  extension  of  the  previous  computations  would  be  to  determine  the  potential 
coefficients  from  5°  mean  anomalies  by  collocation.  However,  as  the  number  of 
observations  then  increases  from  416  to  1654,  it  is  no  longer  a standard  procedure 
to  invert  the  auto-covariance  matrix  of  the  system.  Most  computers  cannot  even 
store  such  a large  matrix.  We  are  going  to  estimate  the  necessary  computer  time 
as  follows.  In  Table  5 the  total  number  of  necessary  multiplications  f (M)  for 
computing  the  potential  coefficients  by  collocation  to  degree  n0  for  M observations 
are  given.  The  direct  method  implies  that  the  method  of  Cholesky  is  used  for  the 
inversion  of  the  auto-covariance  matrix.  Let  us  assume  that  the  total  computer  time 
T (M)  is  proportional  to  f (M).  (We  disregard  the  time  needed  for  addition  operations.) 
Then  we  obtain  from  Tables  4 and  5 for  no  = 20  (no  elevation  information  included): 

T (1654)  = W-6)  f l1654)  = 3*  50* 

' f (416) 

Thus  the  necessary  computer  time  is  so  large  that  we  should  really  consider  whether 
it  could  be  reduced  by  modifying  the  method. 


Method 

Degree  of  Expansion 

CPU  Time 

Integration 

25 

19s 

Collocation 
| with  Elevations 

20 

12m  13s  | 

Collocation 
without  Elevations 
(COVA  in  Table) 

20 

9m  21s 

1 

Computer:  IBM  370/168 
No.  of  Observations:  416 
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One  possibility  might  be  to  determine  the  coefficients  and  their  accuracy 
estimates  according  to: 


,C  n«-. 

hcT 

(6.1) 

is  } 

° na 

- u 

Ag 

and 

— kt 

mc 

c/hc 

(6.2) 

{-2} 

'■ms  ; 

= Co  - 

UsThsJ 

where  the  vectors  of  unknowns  (hc  and  hs)  are  given  by: 


(6.3) 


(C  + D) 


For  each  coefficient  to  be  determined,  the  weights  (hc  or  hs)  can  be  computed  iteratively 
by  the  following  formula  (cf.  Miller,  1974): 


(6.4) 


h"’>  + 0 [c  - <C  + D,  h1"*) 


Ml 


M 1 


M 1 


M 1 


where 


0 < $ < 2/ A 

A„»x  = maximum  eigen  value  of  C + D 
V = iterative  step:  0,  1,  2,  ... 

M = number  of  observations  (mean  anomalies) 

The  starting  value  h^  for  the  iteration  is  most  conveniently  given  in  the  spherical 
approximation.  By  assuming  that  all  mean  anomalies  are  located  on  the  mean  earth 
sphere  of  radius  r,  we  arrive  at  the  following  formula  from  (3.  lb)  after  replacing 
Ag  and  C by  4g  and  /S„  Cna,  respectively  [cf.  Rapp,  1977a,  formula  (30)]: 

M 

(6.5a)  CnB=V  hi0)Aik 

i 

k = l 
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where 


(6.5b) 


1 

4 V £n  ( r s/r  n,)!,+ 3 (n-1) 


(sin  <p)  cos  m A. 


Aak 


d a 


(o)  Substituting  cos  m A under  the  integral  by  sin  m A we  obtain  the  weights 
(hk  ) for  and  coefficients  Su  in  (6.5b)  and  (6.5a),  respectively.  Even  simpler 
approximations  are  obtained  for: 


(6.6) 


r \n+s_  _ eCosmA. 

— ) Pn»(sin«3)i  !•  A<7k/[4ny(n- 1)] 

r ® ' '•sin  mA-1 


which  formula  is  given  from  (6.5b)  by  the  approximation: 

cosmA^  _ _ cosmA 

, } da  ss  pn,  (sin  to)  j - 
sinmA*1  ^sin  mA 

where  <0  and  A are  given  in  (3.  7 ).  As  the  elevation  of  the  highest  mountain  is  less 
than  0.2%  of  the  mean  earth  radius,  we  can  expect  that  the  iteration  error  in  (6.4)  is 
insignificant  after  a few  iterations.  Again,  it  should  be  emphasized  that  the  spherical 
harmonic  coefficients  so  determined  refer  to  the  Bjerhammar  sphere  (of  radius  rB) 
and  should  be  multiplied  by  (rB/a)n+2  in  order  to  be  consistant  with  other  coefficient 
determinations,  which  usually  refer  to  the  sphere  of  radius  a = 6378.140  km. 

In  the  approximate  formula  (6.  5b)  and  (6.  6)  we  have  disregarded  the  noise 
covariance  function  D.  When  considering  the  noise  covariance  function: 


- r 

Act.  J 


Pn 
Aok 


(sin 


V)  { 


d (P,Q)  = Y d „ (rB‘7rp  rQ)"+a  Pn  (cos 

n=a 


the  following  weight  function  can  be  derived  for  the  spherical  case  [see  the  Appendix, 
formula  (A. 9)]: 


to)  t r ^ ^ c c m A . 

(6*7)  hk  = ( — J — -A-  PnB  (sin  <p)  j - [ Aok/4ny(n-l) 

\rs/  cn+dn  *•  s in  m A ^ 


In  the  same  way,  if  we  assume  that  the  errors  between  the  blocks  are  uncorrelated, 
the  following  weight  function  can  be  derived  [formula  (A.  11)]: 
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£n 

rB 

r m 

where 

e 2 = estimated  mean  square  error  of  the  observation  in  block  k. 

By  using  the  iterative  formula  (6.4)  we  avoid  the  inversion  of  the  auto-covariance 
matrix.  Formula  (6.5)  is  theoretically  attractive  in  , because  it  implies  that  the 
iterative  collocation  is  carried  out  with  the  solution  of  the  integration  method  as  the 
original  step. 

Finally,  we  compare  the  number  of  necessary  matrix  • ^rations  for  computing 
the  coefficients  and  their  accuracy  estimates  by  direct  colhx  ation  [formulae  (3.  4)  and 
(4. 3 a-b)]  and  the  proposed  iterative  method  [formulae  (6.1),  (6.2)  and  (6.3)].  In  the 
direct  method  the  computations  of  the  matrix  inverse  and  the  accuracy  estimates  are 
the  most  laborious  operations.  For  the  comparison  we  assume  that  h^  of  (6.4)  is 
a priori  given  and  that  v0  steps  are  necessary  in  the  iterative  method.  The  numbers 
of  necessary  operations  for  a determination  of  the  accuracy  estimates  to  degree  n0 
[i. e.  for  (n0+  l)3  coefficients]  are  summarized  in  Table  5 (the  direct  method  in 
accordance  with  Westlake,  1968,  Table  7.1). 

Table  5 


Number  of  Necessary  Matrix  Operations  to  Compute 
Accuracy  Estimates  to  Degree  n0 


(6.8) 


(°) 


Act, 


ik  - 


y(n-l)4n 


© 


_£jl 


. cos  m a 


3(r.+2) 


^ + (2n+l)  e 2 


Pnm(sinp)t)  j 


L sin  m A - 


Operation 

Direct  Method 

Iterative  Method 

(C  hole  sky) 

[formulae  (6.2)  and  (6.3)] 

Addition 

M3-2M2+  M+  (M3+  M)(n0+ 1)2 

{(M  + l)  Mtb+M}(n0  + 1)2 

Multipl  ication 

£m3+|m2-M+  2M2  (n0+l)2 

(M2F0  + M)(n0+  l)2 

Vo  - number  of  iterative  steps 
M = number  of  mean  anomalies 
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As  the  multiplications  are  the  most  time-consuming  operations,  we  limit  the 
following  comparison  to  those  figures.  Then  we  obtain  from  the  table  that  the 
iterative  method  is  more  efficient  whenever-. 


(M  V0  + 1)  M (n0+  1 f<  a M3  + | Ma  - M + 2M2  (n0  + l)2 


From  this  inequality  we  obtain 


(6.8) 


n0  < 


+ 3M-2 
2M(V0-2)+2 


or,  approximately,  for  V0>  2 


(6.8  ) 


fi 


M+  3 


2(^-2) 


- 1 


Formula  (6.8)  is  illustrated  in  Table  6. 


Table  6 


The  Maximum  Integer  (n0)  Satisfying  (6.8) 
for  Various  M and  V 0 


M V0 

2 

3 

4 

5 

10 

100 

70 

6 

4 

3 

1 i 

416 

294 

13 

9 

7 

4 i 

1000 

707 

21 

14 

11 

6 i 

1654 

1169 

27 

19 

15 

9 1 

5000 

3535 

49 

34 

27 

16  | 

10000 

7071 

69 

49 

39 

24  j 

n0  = maximum  degree  of  series  expansion 
M = number  of  observations 
v0  = number  of  iterative  steps 


The  table  shows  that  the  iterative  method  is  favorable  merely  for  up  to  2 or  possibly 
3 necessary  iterations.  However,  as  earlier  discussed,  the  available  approximations 
h(°)  could  very  well  meet  such  a requirement. 
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7.  Conclusions 


Geopotential  coefficients  determined  by  collocation  were  found  to  agree  some- 
what better  with  the  GEM  9 coefficients  than  the  coefficients  determined  by  pure 
integration  of  10°  mean  anomalies.  This  result  is  probably  due  to  the  incorporation 
of  a weighting  of  the  observations  with  respect  to  their  a priori  accuracies.  However, 
this  gain  is  achieved  at  the  cost  of  several  times  more  computation  time. 

By  the  inclusion  of  the  elevation  information  in  collocation,  the  coefficients 
to  degree  20  changed  by  3%  on  the  average.  A surprisingly  large  difference  of  10% 
was  obtained  for  n = 3.  The  RMS  changes  of  the  undulation  and  the  anomaly  were 
2 meters  and  0. 7 mgal,  respectively.  In  general,  however,  we  might  expect  that  the 
10  blocks  give  a too  rough  approximation  to  the  topographv  to  re-  eal  any  more 
significant  magnitudes  of  the  changes  of  the  coefficients  when  . luding  a correction 
for  the  topography  (the  Molodensky  G!  term).  A possible  explanation  of  the  10% 
differences  for  n = 3 in  various  methods  might  be  the  non-symmetric  distribution  of 
the  continents  between  the  northern  and  the  southern  hemisphere. 

From  the  comparison  of  the  coefficient  accuracy  estimates  for  various 
methods  (Table  3)  we  conclude  that  the  error  covariances  between  the  mean  anomalies 
should  be  taken  into  account  in  the  computations  by  direct  integration.  In  collocation 
it  seems  less  important  to  include  these  covariances  in  the  computations. 

0 ^ natural  continuation  of  the  above  study  would  be  to  compute  the  coefficients 

for  5 mean  anomalies  by  collocation.  However,  due  to  the  difficult  task  to  invert  an 
auto-covariance  matrix  for  more  than  1600  observations,  the  original  method  should 
irst  be  modified  according  to  the  iterative  method  described  in  section  6.  The  method 
includes  the  solution  by  integration  as  a preliminary  step.  As  this  technique  avoids 
the  inversion  of  the  auto-covariance  matrix,  a considerable  gain  in  computer  time  can 

be  expected.  Another  possibility  would  be  to  determine  the  auto -covariance  matrix 
in  an  iterative  way. 
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Appendix 

It  is  desired  to  determine  the  spherical  harmonic  coefficients 
(referring  to  the  Bjerhammar  sphere  of  radius  rB)  implied  by  a mean 
anomaly  field  Ag  on  a sphere  of  radius  r,.  The  following  covariance 
functions  are  given: 


(A.l) 


ce  \ (p)  = / cov  (CnV.  2i) 
c.  J ! cov  (S„V, 


Y,'  »'(P) 


where 


(A.  2) 
and 
(A.  3) 


W = cn^/ /y(2n’ + l)(n’-l) 

v.v(P)  - 

00 

c(Q,P)  = cov(2S&,5&)  =Yc„^(rBVr9rpr-Pn(co8^p) 


n =a 


k 

d(Q,P)  = cov  ( Cg  , eP ) =Y  dnfln2(r03 /rQ  rPy>+a  Pn(cos(pqp) 


The  solutions  for  Cnvand"5rVby  the  method  of  least  squares 
collocation  are  given  by  (see  Sjoberg,  1978) 


(A.  4) 


_l  f r ;hc  <Q) 

4r  J J t hs  (Q) 


} 4g  (Q)  dc^ 


where  the  weight  functions  he  (Q)  and  hs  (Q)  are  given  by  h(Q)  of  the  following 
Wiener- Hopf  integral  equations 


(A.  5)  ill  } (P)  = JJ  ^Q)  {C  (Q’  P)  + d(Q’  P)  } da« 

A A 

Proposition!:  The  weight  functions  for  C„.  and  Snl  in  (A. 4)  implied  by 
(A.  1)  - (A.  3)  and  (A.  5)  for  n>  = rQ  = r,  = constant  are  given  by: 

(A.6)  h(Q)  = + YB1(Q)/y/3n(n-l) 

\ rB  / cn  + dn 
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I -«> 


Proo^  We  expand  (A.  2)  and  (A.  3)  accordingly: 


? 

I I 


(A.  2') 

c(q,p)  »y  y 

a =8 

/ 

2n+ 1 ' 

' re  \2("  +a  ) 

,~tj  Yn.(Q)  Yal(P) 

and 

(A.  3’) 

«*(Q.p)  - y y 

■ — * 

n=3  ■*- 

dAa 

n 2n  +1 

(f-)a(n+a)  Yn" (Q)  Yn.<p) 

(A.  7) 

h JJ  Y-  Y»  - 

f 1 if  n = p and  m = q 
\ 0 otherwise 

Inserting  (A.l),  (A. 2'),  (A. 3')  and  the  expansion 


h<y>  = t t hn*  yb,(q)  i 

n=3  a=-n 


Corollary  1;  For  d(Q,P)  of  Proposition  1 replaced  by 
d(Q,P)  = e3(Q)  6(iA,p) 
where  6(^p)  is  the  Dirac  delta  function,  we  obtain 


(A.  10)  h(Q)  = _ aftn  (l g_  \ 

y(n-l)  \r„  / 


Nn+3  J 

cn/9nVB/rB)3'n+3)  + (2n+l)  es(Q)  Ynit(Q) 


Proof:  As 


^ JJ  f3  (Q)  6(^,P)  Yna(Q)  doQ  = ea(P)Ym(P) 
formula  (A.  8)  becomes  in  this  case 

')n'(rB/r/+J  Yn'»'(P)  = Yhna  J8n2  (^)a(n+a)+e3(P))  Yn,(P) 

n>m 

and  the  proof  follows  after  noting  that 

h na  = 0 for  n ^ n ' and  m / m! 

In  accordance  with  formula  (A. 9)  we  obtain  in  this  case: 


(A.  11) 


j_  rr 


II  h<Q> 


dcr 


/'JkY+a 1 v 'Aok 

4n  JJ  y(n-l)  ' r,/  c n /9n3  ( rB/r„ )3 ("  +s ) + (2n+l) e^Qk ) ^ n* (Qy^4rr 


Proposition  2;  The  error  estimates  for  Cn„  and  SnB  of  (A.  4)  and  (A.  5)  are 
given  by 


= Co_Ii>  JJIJh(Q)  b(Q)  {c(Q,Q')+d(Q,Q')  } dcTqd<v 

where  C0  is  the  a priori  variance  of  the  coefficients. 

Pro°ft  We  consider  only  the  estimate  Cn„. 

= E f(Cn.-CnB)3}  = E {C„3 }+  E fCni  } - 2 E [CMCnB  } 
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where 


E { CB.a  } = Co 

E { C„a  } = ^-5  JJ  j j h(Q)  h(Q ')  E fTSfe,  ,}  do,  do,  ' 

= leg'll  J J^Q)  h(Q'>  {c(Q,Q')  + d(Q,Q')  }do,do5}, 
E {Q.Q.  } - ~ J j*  h(Q)  E f Cni  2^  } da, 


= ^rIIh(Q)  Cc(Q)  dC* 

= lfrT  IJ  I Ih<^>  h<Q)  (C(Q,  Q ) +d(Q,  Q)  } da^,  do,. 
From  these  deductions  the  proposition  follows  immediately. 
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Table  A.  1:  Fully  Normalized  Geopotential  Coefficients  Determined 
from  416  10°  Mean  Anomalies  by  Collocation.  No  Mean  Elevations 
Included.  All  Coefficients  multiplied  by  106.  a = 6378140  m. 


n 

n 

C 

B 

81CTA 

BICTtA 

F 

n 

C 

B 

BICttA 

SfCffA 

2 

-403.4449 

0.2016  . 

2 

i 

9.3000 

-0.0996 

0. 1790 

0.  170B 

2 

2 

2.4400 

-1.3057 

0. 1670 

0.  1597 

3 

9.7 133 

0.2304 

n 

i 

1 . 0349 

0. 1234 

0.2263 

0.2261 

9 

2 

1.1793 

-9.4719 

0.2263 

0.2243 

• 

3 

0.7301 

I.G3I0 

0.2193 

0.2203 

4 

0.0957 

0. 1324 

4 

t 

-9.4544 

-0.4057 

9. 1219 

0. 1217 

4 

2 

0.3170 

9.3939 

9. 1223 

9. 1227 

4 

3 

9.7935 

-0.3966 

0. 1 160 

0. 1 100 

4 

4 

-0.2268 

9.4  66 

0.1151 

9. 1141 

3 

0.2234 

0.0096 

s 

1 

-0. 1576 

-0.2402 

0.0021 

9.0017 

0 

2 

0.4107 

-0.2016 

0.0020 

9.0024 

s 

3 

-9.2440 

-0. 1369 

9.001  1 

o.oau 

8 

4 

-9.9921 

-9.0276 

9.0771 

9.0769 

0 

5 

9. 1064 

-0.4974 

0.0740 

0.0766 

6 

-0. 1376 

9.0340 

4 

1 

9.1120 

0.9633 

9.0466 

0.0437 

4 

2 

9.2019 

-9.3649 

9. 0464 

9.0459 

4 

3 

-0.0173 

-0.0623 

0.0439 

0.0433 

6 

4 

-0. IUOO 

-9.4041 

0.0426 

9.0423 

4 

0 

-0.3769 

-0.3370 

9.0363 

0.0396 

4 

6 

9.9997 

-9.2402 

9.0371 

0.037G 

7 

0.2002 

0.0347 

7 

1 

0.2010 

0.0466 

9.0313 

9.0310 

T 

2 

9.21150 

t.  1194 

9.0497 

9.0499 

7 

3 

9. 1615 

-0. 1003 

9.0301 

0.0301 

7 

4 

-9. 1372 

-0. 1660 

9.9490 

0.9409 

7 

G 

9.0096 

-0.0060 

0.9466 

0.0472 

7 

6 

-9.3070 

0.2000 

9. 0440 

0.9449 

7 

7 

-9.0323 

-0.0921 

9.0459 

9.0449 

0 

9.0422 

9.0441 

u 

1 

-9.0423 

0.0373 

0.0402 

0.0409 

n 

2 

9. 1333 

0.0993 

9.9399 

9.9404 

n 

3 

9.0704 

-0.0234 

9.0300 

9.0303 

B 

4 

-9.2167 

9.0313 

9.0392 

9.0307 

u 

3 

-0.9174 

0.0349 

0.0366 

0.0373 

0 

6 

-9.9006 

0.  1900 

0.0338 

0.9334 

0 

7 

0.0303 

0.0713 

9.0330 

0.0336 

0 

(1 

-9. 1417 

9.0911 

9.0343 

9.0330 

4 

9. 1336 

0.0409 

9 

1 

9. 1723 

-0.0220 

9.0302 

9.0301 

9 

2 

9.1171 

-9.9043 

9.03110 

9.0304 

4 

3 

-9. 1706 

-0.0303 

0.0373 

0.0372 

9 

4 

-9.039(1 

0.9490 

0.0360 

9.0363 

9 

0 

-4.0572 

-0.0336 

0.0362 

9.0367 

9 

4 

9.01114 

0. 1023 

9.0347 

0.0344 

9 

7 

-9.01123 

-0.4210 

9.0336 

0.0336 

9 

0 

9.2930 

9.0063 

9.0327 

0.9324 

V 

9 

-0.9159 

0.0237 

0.0333 

0 . 0333 

10 

9.9303 

0.0331 

10 

1 

9. 1146 

-0.0147 

0 . 0326 

0.932(1 

19 

2 

-9.0637 

-9. 0487 

9.0324 

9.9326 

l« 

3 

-0.9416 

-0.0954 

0.4310 

9.0319 

19 

4 

-9.0674 

-0.9963 

9.0309 

9.9303 

10 

B 

-0.0261 

-0.9113 

0.0301 

0.9306 

to 

6 

-9.047G 

-0.0479 

0.0299 

0.9296 

10 

7 

9.0024 

-0. 0 1117 

0.0270 

9.9274 

19 

II 

9.0413 

-9.9333 

0.0273 

0.0270 

10 

9 

0. 1013 

-0.0222 

0.0263 

9.9263 

19 

19 

9. 1159 

-9.9309 

9.9277 

9.0274 

II 

-0.9933 

9.0336 

II 

1 

-9.9214 

0.0107 

0.0310 

9.0311 

II 

2 

-9.9219 

-9.0922 

9.9317 

9.0310 

II 

3 

-0.0747 

-0. 1 126 

9.0306 

0.0304 

II 

4 

-0. 1044 

-9.0707 

9.0302 

9.0299 

II 

0 

0.0123 

0. 0074 

0.02IU1 

0.0292 

11 

9 

-0.9112 

0.0100 

0.0209 

0.0207 

II 

7 

0.0373 

-0. 1030 

0.0279 

0.02110 

II 

n 

0.0992 

9.0613 

9.0266 

0.0263 

II 

9 

-0.9456 

9.0017 

0.0263 

0.0264 

II 

10 

-9.0269 

9.9139 

0.9263 

0.0233 

II 

II 

0.0769 

-0.0164 

0.0271 

0.0271 

12 

9.0330 

0.0321 

12 

f 

-0.0213 

-9.0334 

0.0297 

0.0298 

12 

2 

9.9030 

-9.9399 

9.0303 

9.0306 

12 

3 

0.9097 

9.0307 

0.0293 

0.0294 

12 

4 

-9.0653 

-9.0343 

9.0200 

0.02117 

U 

G 

0.0376 

-0.9023 

9.0201 

9. 0283 

12 

6 

0.0233 

9.0496 

9.9273 

0.0272 

12 

7 

-0.0249 

9.0267 

0.0272 

0.0271 

12 

n 

0.0190 

0.0264 

0.0263 

0.0261  • 

12 

9 

0.9106 

9.0169 

9.0232 

9.9231 

12 

10 

-0.0094 

-9.0263 

9.0234 

9.9233 

12 

II 

9.9009 

9.0064 

0.0230 

9.0231 

12 

12 

0.9134 

-9.0094 

0.0268 

9.0264 

12 

0.0473 

0.0303 

12 

1 

-9.0016 

9.9106 

0. 0200 

9.0292 

19 

2 

9.0003 

-9.9414 

9.9209 

9.9293 

12 

0 

-0.0144 

9.0421 

9.02117 

9.02(16 

13 

4 

9.0919 

-9.9060 

0.9277 

9.0273 

13 

9 

0.9306 

9.0397 

0.4271 

0.0273 

13 

6 

-9.0246 

-9.00(16 

9.0263 

9.0265 

13 

7 

-0.0171 

9.0157 

0.0239 

9.0239 

13 

0 

-9.0237 

9.9200 

0.9230 

9.9253 

13 

9 

-0.9140 

9.0334 

0.0210 

0.0247 

13 

10 

9.0339 

-9.0149 

9.0243 

0.0249 

13 

II 

0.0043 

9.0201 

9.0243 

0.0246 

13 

12 

-9.9016 

9.9943 

9.0242 

0.0248 

13 

13 

-0.0349 

9.0637 

0.0237 

0.0260 

14  9.0000  9.02110 
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Table  A.  1 (continued) 


14 

1 

• .0204 

0.0904 

0.0276 

0.0209 

t4 

2 

-0.9307 

-0.0173 

9.6272 

0.9772 

14 

9 

0.0096 

0.0002 

0.0272 

0.9270 

14 

4 

-0.0916 

-0.0209 

0.0263 

0.0260 

C 

0.0147 

-0.0070 

0.0234 

0.02S4 

14 

6 

-0.09113 

0.00)0 

0.0231 

8.0230 

14 

7 
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0.0241 

14 

n 
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4 
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14 

10 

0.0160 

-0. 00*16 
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14 

11 

0.0336 
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0.0223 

14 
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0 023' 

14 

13 
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14 

16 
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0.01  17 

9.0244 

0.0246 

13 
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0.0266 

IS 

1 
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18 

2 

0.0043 

0.9133 

0.0266 

0.0267 

r 

3 
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0 0452 

0.0262 
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13 

6 

6.0042 

-9  017) 

0.023(1 

13 

6 
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0.0247 
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IS 
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- QtJ  * 
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13 
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0.0303 
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18 
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9 . 023  1 

IS 

9 
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IS 
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Table  A.2:  Fully  Normalized  Geopotential  Coefficients  Determined 
from  416  10°  Mean  Anomalies  by  Collocation.  Approximate  Mean 
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Table  A.  3:  Fully  Normalized  Geopotential  Coefficients  Determined 
from  416  10°  Mean  Anomalies  by  Collocation.  No  Elevations  Included. 
Noise  Covariance  Matrix  D = 0.  a = 6378140  m.  Coefficients 
multiplied  by  106. 
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Table  A.  3 (continued) 
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Table  A.  4;  Fully  Normalized  Geopotential  Coefficients  Determined 
from  416  10°  Mean  Anomalies  by  Integration.  All  coefficients  mul- 
tiplied by  10s.  (Error  Estimates  not  corrected  by  AT1.) 
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